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ABSTRACT 


Alkaline protease producing microorganism's isolated total of 74 and 69 fungi from beach soils and black soils respectively were 10 from Kakinada beach soil and 9 
from Vishakhapatnam beach soil and 40 isolates from Tenali black soil. The best producing organism identified and characterized as Pencilium oxalicum KRSS-S- 
F P10. Taguchi was used to know active and independent variables screened by using PDA media variables was employed to derive statistical model (Response surface 
methodology: RSM) for optimizing the medium composition. The alkaline protease activity was in between 2.0U/mL to 10.76U/mL and The Delta value indicates the 
effect of component in the form of ranking has given to the variables. To confirm these results, experiments were carried out using RSM these nutrients with same 
concentrations and it was observed that the mean value of alkaline protease was 10.13U/mL as compared to the predicted value of 10.76U/mL using MINITAB. The 
final optimized medium has given approximately three-fold increase in alkaline protease production in comparison with one-factor-at-a-time (5.12U/mL) method. 
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INTRODUCTION: 

Currently enzymes can be used in basic and applied arenas of research as well as 
ina wide range of product design and manufacturing processes, such as food, bev- 
erage, pharmaceutical, detergent, leather processing and peptide synthesis 
industries with the estimated value of the global sales over 3 billion USD (Jordan 
Chapman ef al., 2018). Of the whole industrial enzymes, 75% are hydrolytic. 


Three soil samples were collected in depth of 5-6 cm from soil. Alkaline pH was 
observed across the sampling sites and a total of 3 samples were found with 
plenty of micro flora occurrence (Ashok Pandey et al., (2013)). The total micro- 
organism counts of the soils were estimated by standard dilution plate technique. 
The isolated microbes were identified by their cultural and fungi and 
actenomycetes were obtained from the positive three soil samples (Dilution 
Worksheet and Problems (2021)). These are 10 fungal forms from Kakinada, 9 
from Vishakhapatnam and 40 from Tenali black field morphological characteris- 
tics. Total 74 microorganism forms (includes bacteria soils. Alkaline proteases 
fungi were isolated using milk agar plate assay consists of 0.5% casein from dif- 
ferent soils collected from Kakinada, Vishakhapatnam beach soil and Tenali 
black soil fields of Andhra Pradesh. Three soil fungal isolates were examined for 
protease producing microorganism. Isolation and RAPD-identification of the 
highest alkaline protease producer under submerged fermentation, using PDA 
media. 


MATERIALS AND METHODS: 

Reagents: Unless otherwise stated, the chemicals and medium ingredient used 
in this study were purchased from Sigma Chemical Co. (St. Louis, Mo, USA) and 
BD Bioscience (Le Pont de Claix, France). All other chemicals were also of ana- 
lytical grade. 


Optimization of enzyme production: Characterization of the different factors 
for alkaline protease production was optimized by applying RSM. Then statisti- 
cal model was obtained using Central composite design (CCD) with three inde- 
pendent variables alkaline protease (X,), (X,) and MnSO, concentration (X,). 
CCD maximizes the amount of information that can be obtained while limiting 
the number of individual experiments (Kunamneni and Singh, 2005). Each fac- 
tor in the design was studied at five different levels (Table-1). A set of 20 experi- 
ments were performed. All variables taken at a central coded values were listed in 
Table-2. Upon completion of experiments, the average of alkaline protease pro- 
duction was taken as the dependent variable or response. 


pH Casein Glucose FeSO, NaMoO, MnCl, 


— 
— 


Table 4.5.1.1. Optimization of the factors of Taguchi L,, 
CaCl, 


Statistical analysis: One factor ANOVA in Taguchi and MINITAB software was 
used to analyse the experimental results. Statistical analysis is considered as an 
important tool to interpret and summarise the experimental data. All the experi- 
ments results were carried out in duplicate and data obtained is represented as 
mean + SD in this paper. P-values of the experimental data are reported respec- 
tive figures. The P-value less than 0.05 is considered to be statistically signifi- 
cant. 


RESULTS AND DISCUSSION: 

4.5.1 Optimization of the factors of Taguchi L,,: 

Using L,, orthogonal array design approach, the relationships between medium 
component variables and their concentrations could be worked out. The optimal 
combinations and the concentration of the factors required to achieve the highest 
alkaline protease activity was represented in Table 4.5.1.1. The alkaline protease 
activity was in between 2.0 U/mL to 10.76 U/mL at different levels of various 
independent factors. The Taguchi approach suggests to analyze variation by 
using an appropriately chosen signal-to-noise (S=N) ratio and also serve to ana- 
lyze the average response for each run in the inner array. The Table-4.5.1.2 and 
4.5.1.3 represent the response table for mean and S/N ratio to understand the 
delta and the rank value of the system. The Delta value indicates the effect of that 
component whereas rank (Table-4.5.1.4) based on the delta values serialized the 
factors from the greatest effect to the least effect on the response. 


In the present study, it can be seen that for each of the seven variables at three lev- 
els, one level increases the mean compared to the other level (Table-4.5.1.5). 
Thus the factors NaNO, at level-1, MgSO, at level-1, ZnSO, at level-1, TA at 
level-1, KH,PO, at level-2, FeSO, at level-1, KCI at level-1 shows a main effect. 
These levels also represent the optimal concentrations of the individual compo- 
nents in the medium. The Fig.10 and Fig.11 represent the main effect plot of 
means and S/N ratio in the system. Final optimized levels of each factor in the 
medium for alkaline protease production had shown in Table-4.5.1.1.To confirm 
these results, experiments were carried out using these nutrients with same con- 
centrations and it was observed that the mean value of alkaline protease was 
10.13 U/mL as compared to the predicted value of 10.76 U/mL using MINITAB. 
The final optimized medium (Table-4.5.1.6) has given approximately three-fold 
increase in alkaline protease production in comparison with one-factor-at-a-time 
(5.12U/mL) method. Therefore, it should be considered that the selected condi- 
tions were the most suitable in practice and NaNO,, casein and KCl were identi- 
fied as the major influencing factors in the enzyme production. 


Mgso, | NH,n0, KH,PO, 


rt 


— 
N 


1 
2) 
3 
4 
5) 
6 
7 
8 
9 


= |S) oS fs || OS | Bs 
elalalalelajela 


DIADININIAI AIAN] 


elafalalelalela 
NIN Nl aAlaAla|y| ry 


International Education & Research Journal [IERJ] 


Research Paper 


E-ISSN No : 2454-9916 | Volume: 8 | Issue: 6 | Jun 2022 


Glucose 
40.6685 


Casein 
43.5218 


Level 


40.4238 25.1055 


Table 4.5.1.2. Response Table for Signal to Noise Ratios 
NaMoO, 
28.9432 


MnCl, 
22.9226 


CaCl, 
18.0618 
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Fig. 11: Main Effects Plot for Means 
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Fig. 12: Main Effects Plot for SN ratios 


4.5.2. Optimization of the selected medium components by RSM: 

The Taguchi results showed that among seven different culture components, four 
have significant impact on enzyme production (Mishra, Rashmi, (2020)). 
Response surface methodology was used to investigate the effects of four com- 
ponents viz. MnClL,, Casein, NaMoO, and pH by taking as independent variables 
onalkaline protease enzyme production (Malhotra and Chapadgaonkar, 2020). It 
was found that alkaline protease production was enhanced up to 10.76 U/mL in 
the medium optimized by RSM. 


The CCD matrix in terms of code and actual values of independent variables is 
given in Table-4.5.2.1. The second order model equation was quantified by using 
Analysis of variance (ANOVA). ANOVA of the quadratic model indicated that 
the model is significant. In this study it was found that A, C, D, A2, H2, AB, AC, 
AH, BH and CH (A=MnCl, 0.005 g/L, B=casein 13.1521g/L, C= NaMoO, 
0.0205 g/L, and D=pH 9.0/10.0) are the model terms. The second order response 
model was made after analysis of regression. The model can be illustrated as fol- 
lows. 


Final Equation in Terms of Coded Factors: 


R1=+523.9752.65*A3642.48*B45.32*C173.27*D+79.31*AB183.06*AC1 18. 
23*AD+77.81*BC+842.31*BD101.88*CD+123.35*A‘%2+2398.65*B%2+1 17. 


8) 


52*C42+200.31*ABD158.81*ACD+181.31*BCD130.35*A%2D396.40*B%2 
D138.66*C’2D. 


Final Equation in Terms of Actual Factors (pH 9.0): 
R1=+1.22166E+0051.83737E+006*MnCl21884.72947*Casien+ 14828.82972 
*NaMo042420.00000*MnCl2*Casien2.487 18E+005*MnCl2*NaMo04530.7 
6923*Casien*NaMo04+ 1.01484E+007*MnCl2“2+27.9505 1 *Casien*2+6.73 
725E+005* NaMoO4“2 and 


Final Equation in Terms of Actual Factors: pH10.0: 
R1=+37550.07733-76878.05410*MnC12-2039.89282*Casient+2.87974E+ 
005*NaMo004+5592.50000*MnCl12*Casien3.50641E+006*MnCl2*NaMo04 
+1328.84615*Casien* NaMo042.79959E+005*MnC12%2+20.02245*Casien* 
25597.90761*NaMoO04“2 


Table 4.5.2.1: Actual and predicted values 


Internally Externally Influence on 


Run Actual Predicted Studentized Studentized Cook's Fitted Value Standard 
Order Value Value Residual Leverage Residual Residual Distance DFFITS Order 
1 8808.00 8484.71 323.29 0.670 0.723 0.714 0.053 1.018 5 
2 702.00 697.24 4.76 0.166 0.007 0.007 0.000 0.003 15 
3 709.00 1525.14 -816.14 0.607 -1.675 -1.760 0.217 -2.189 * 10 
4 350.00 350.70 -0.70 0.166 -0.001 -0.001 0.000 -0.000 39 
5 386.00 350.70 35.30 0.166 0.050 0.048 0.000 0.022 37 
6 370.00 618.29 -248.29 0.607 -0.509 -0.500 0.020 -0.622 29 
i 15850.00 16145.34 -295.34 0.607 -0.606 -0.596 0.028 -0.741 ll 
8 55.00 596.58 651.58 0.670 1.458 1.503 0.216 2.141 * 28 
9 356.00 350.70 5.30 0.166 0.007 0.007 0.000 0.003 36 
10 759.00 697.24 61.76 0.166 0.087 0.085 0.000 0.038 19 
ll 12018.00 10723.21 1294.79 0.607 2.657 3.219 0.546 4.003 * 31 
12 12.00 1060.34 -1048.34 0.607 2.151 -2.391 0.358 -2.974* 12 
13 45.00 -136.68 181.68 0.670 0.407 0.398 0.017 0.567 24 
14 2822.00 3926.26 -1104.26 0.670 2.471 2.889 0.619 4.115 * 26 
15 5308.00 5511.03 -203.03 0.670 0.454 0.445 0.021 -0.634 25 
16 152.00 -609.22 761.22 0.670 1.703 1.795 0.294 2.557 * 8 
17 304.00 538.46 -234.46 0.607 0.481 -0.472 0.018 0.587 33 
18 5208.00 5422.66 -214.66 0.670 0.480 -0.471 0.023 -0.671 22 
19 312.00 350.70 38.70 0.166 -0.054 -0.053 0.000 -0.024 40 
20 8305.00 8116.10 188.90 0.670 0.423 0.414 0.018 0.589 1 
21 308.00 43.34 264.66 0.607 0.543 0.533 0.023 0.663 34 
22 45,00 -1037.91 1082.91 0.670 2.423 2.810 0.595 4.002 * 23 
23 701.00 697.24 3.76 0.166 0.005 0.005 0.000 0.002 17 
24 9208.00 8809.37 398.63 0.670 0.892 0.887 0.081 1.263 6 
25 740.00 1516.95 -7716.95 0.607 1.594 -1.663 0.197 -2.068 * 14 
26 332.00 350.70 -18.70 0.166 0.026 -0.026 0.000 -0.011 38 
27 160.00 ~466.83 626.83 0.670 1.403 1.440 0.199 2.050 * 4 
28 63.00 -130.31 193.31 0.670 0.433 0.424 0.019 0.603 27 
29 363.00 350.70 12.30 0.166 0.017 0.017 0.000 0.008 35 
30 708.00 697.24 10.76 0.166 0.015 0.015 0.000 0.007 18 
31 4967.00 5639.93 ~672.93 0.670 -1.506 1.559 0.230 -2.220* 21 
32 8894.00 8537.76 356.24 0.670 0.797 0.790 0.064 1.124 2 
33 40.00 1304.59 -1264.59 0.607 2.595 3.105 0.521 -3.862 * 32 
34 789.00 697.24 91.76 0.166 0.129 0.126 0.000 0.056 20 
35 144.00 ~449.89 593.89 0.670 1.329 1.356 0.179 1.932 * 7 
36 760.00 1326.73 -566.73 0.607 -1.163 -1.174 0.105 -1.460 * 13 
37 147.00 -404.50 551.50 0.670 1.234 1.251 0.154 1.782 * 3 
38 755.00 697.24 57.76 0.166 0.081 0.079 0.000 0.035 16 
39 322.00 43.51 278.49 0.607 0.571 0.562 0.025 0.698 30 
40 777.00 1304.55 -527.55 0.607 -1.082 -1.087 0.091 -1.352 9 


The significance of second order model equation was verified by F (ANOVA) 
test (Table-4.5.2.3.). For this model, if the F test is significant at 5% level i.e. 
P<0.0001, then the model 5 is considered as fit. ANOVA analysis can explain 
whether the more complex model is required for a better fit and it gives the values 
of the model (Rouaa Daou et al., 2021). According to the ANOVO of the qua- 
dratic regression model, it was identified that the present model is highly signifi- 
cant and it is justified by the F test value of 49.79 with low probability value 
(Pmodel> F is 0.0001) (Table-4.5.2.3). The "Pred R-Squared" of 0.9793 is not as 
close to the "Adj R-Squared" of 0.9596as one might normally expect (Abiola 
Ezekiel Taiwo et al., 2020). This may indicate a large block effect or a possible 
problem with the model and/or data. Things to consider are model reduction, 
response tranformation, outliers, etc "Adeq Precision" measures the signal to 
noise ratio. A ratio greater than 4 is desirable (Sunil Chamoli, 2015). A ratio of 
31.246 indicates an adequate signal (Table-4.5.2.4.). The three dimensional plots 
of the responses were drawn to predict the alkaline protease for different test vari- 
able values to understand the interaction among the independent factors (Ram 
Kumar ef al., 2018). The plots of A’, B’, C’, AC, BC, BD were constructed by tak- 
ing response (alkaline protease production) at z-axis against two variables (Fig. 
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12). It was observed that medium containing MnCl, 0.005 g/L, casein 
13.1521g/L, NaMoO, 0.0205 g/L, and pH 9.0/10.0. The predicted results com- 
paratively show 0.5-fold increase with Taguchi (10.76 to 15.82 U/mL) alkaline 
protease activity (Hammami, Bayoudh, and Abdelhedi, 2018). Similarly, the 
effects of parameters such as, relative humidity, pH of the liquid medium, and vol- 
ume of inoculums were positive and smaller than the error limits and the change 
in mean effect was also small, in which case according to the decision-making 
procedure of the evolutionary operation technique it was advisable to select a 
new search region and start a new phase of experiments reported (Dhayalan, 
Velramar, Govindasamy, 2022). RSM has been applied for designing of experi- 
ments to evaluate the interactive effects through a full 31 factorial design and 
reported the similar results by Gilmour, Steven, in the year (2006). The optimum 
conditions were casein concentration, 3.22%; fermentation period, 96h; temper- 
ature, 30°C; and pH 9.0. The high value of the regression coefficient (R'=0.9793) 
indicate excellent evaluation of experimental data by second-order polynomial 
regression model. The RSM revealed that a maximum alkaline protease produc- 
tion of 15.82 U/mL was obtained at the optimum conditions and found the similar 
results with Gomaa, (2013). 


Table 4.5.2.2: ANOVA for Response Surface Reduced cubic model 
Analysis of variance table. 


ANOVA for Response Surface Reduced Cubic model 
Analysis of variance table [Partial sum of squares - Type II] 
F 


Sum of Mean p-value 
Source Squares df Square Value Prob > F 
Model 5.722E+008 19 3.012E+007 49.79 < 0.0001 significant 
A-MnCl2 75716.45 1 75716.45 0.13 0.7272 
B-Casien 3.624E+008 1 3.624E+008 599.15 = 0.0001 
C-NaMoO4 = 56104.35 1 56104.35 0.093 0.7638 
D-pH 3.610E+005 1 3.610E+005 0.60 0.4488 
AB 1.006E+005_ 1 1.006E+005 0.17 0.6877 
AC 5.362E+005 1 5.362E+005 0.89 0.3577 
AD 3.818E+005 1 3.818E+005 0.63 0.4362 
BC 96876.56 1 96876.56 0.16 0.6932 
BD 1.938E+007 1 1.938E+007 32.04 = 0.0001 
cD 2.835E+005 1 2.835E+005 0.47 0.5014 
A“2 4.386E+005 1 4.386E+005 0.73 0.4046 
B“2 1.658E+008 1 1.658E+008 274.17 <= 0.0001 
c*2 3.981E+005 1 3.981E+005 0.66 0.4268 
ABD 6.420E+005 1 6.420E+005 1.06 0.3152 
ACD 4.035E+005 1 4.035E+005 0.67 0.4237 
BCD 5.260E+005 1 5.260E+005 0.87 0.3622 
A‘“%2D 4.898E+005 1 4.898E+005 0.81 0.3789 
B“‘2D 4.529E+006 1 4.529E+006 7.49 0.0127 
c*2D 5.542E+005_ 1 5.542E+005 0.92 0.3499 
Residual 1.210E+007 20 6.048E+005 
Lack of Fit 1.209E+007 10 1.209E+006 1192.90 <0.0001 _ significant 
Pure Error 10132.17 10 1013.22 
Cor Total 5.843E+008 39 


Table 4.5.2.3. Response curve standards with factor 1 and R1 


I cl cl 
nt 

[Factor [Estimate |df [Error [Low [High | VIF | 
[Intercept_|523.97__ [1224.29 [56.12 [991.82 |__| 


A-MnCl2 | -52.65 148.81 -363.06 257.76 1.00 
B-Casien | -3642.48 148.81 -3952.89 |-3332.07| 1.00 
Cc- -45.32 148.81 -355.73 265.09 1.00 
NaMoO4 


[D-pH___|-173.27__|1__—*224.29 __|-641.13__|294.58 13.33 _| 


[AB (79.31 [1194.43 | -326.26 1484.88 [1.00 _| 
[AC | -183.06 |1 1194.43 | -588.63 1222.51 [1.00 _| 
[AD -118.23, [1 [148.81 |-428.64 [192.18 [1.00 _| 
IBC (77.81 | 194.43 | -327.76 1483.38 [1.00 _| 
IBD |842.31_|1 [148.81 |531.90 __]1152.73 [1.00 | 
ICDs -101.88 |i 148.81 [-412.29 1208.53 [1.00 _| 
A‘“2 123.35 144.86 -178.82 425.53 1.02 
lo*2_ 7.52 fn [144.86 |-184.65 [419.70 [1.02 | 
Rep —[aseai [tists —[ sens [ate re Lo 
ACD -158.81 194.43 -564.38 246.76 1.00 
IBCD [181.31 |i 194.43 | -224.26 [586.88 [1.00 _| 
B“°2D -396.40 144.86 -698.58 -94.23 1.67 
1 


C*2D -138.66 144.86 -440.84 163.51 1.67 
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CONCLUSION: 

Present investigation aims to develop new bioprocess and economically feasible 
media by using optimization of production medium composition. It effects the 
cost of production medium and found effective production of enzyme. By taking 
PDA media and serially optimized the media by using some statistical tools are 
Taguchi, RSM and Minitab reduces the unwanted variables in the media for good 
fermentation, this step reduces the cost of production. As a result, use of depend- 
ent and independent variables enhance the production of alkaline protease shake 
flask fermentation. The bioprocess optimization predicted results comparatively 
show 0.5-fold increase with Taguchi to RSM (10.76 to 15.82 U/mL) alkaline pro- 
tease activity and also have prepared different production media for enzyme fer- 


International Education & Research Journal [IERJ] 


Research Paper 


mentation. 


REFERENCES: 


18 


XH. 


XIV. 


Dhayalan, A., Velramar, B., Govindasamy, B. (2022). Isolation of a bacterial strain 
from the gut of the fish, Systomus sarana, identification of the isolated strain, opti- 
mized production of its protease, the enzyme purification, and partial structural char- 
acterization. J Genet Eng Biotechnol 20, 24. 


Gilmour, Steven. (2006). Response Surface Designs for Experiments in 
Bioprocessing. Biometrics. 62. 323-31. 


Gomaa E Z. (2013). Optimization and characterization of alkaline protease and 
carboxymethyl-cellulase produced by Bacillus pumillus grown on Ficus nitida 
wastes. Braz J Microbiol. 2013;44(2):529-537. 


Ram Kumar A, Siva Kumar N, Gujarathi AM, Victor R. (2018). Production of 
thermotolerant, detergent stable alkaline protease using the gut waste of Sardinella 


longiceps as a substrate: Optimization and characterization. Sci Rep. 
2018;8(1):12442. 


Parameswaran Binod, Piyush Palkhiwala, Raghavendra Gaikaiwari, K Madhavan 
Nampoothiri, Arvind Duggal, Kakali Dey and Ashok Pandey (2013). Industrial 
Enzymes - Present status and future perspectives for India; Journal of Scientific & 
Industrial Research; Vol. 72: 271-286. 


Dilution Worksheet and Problems (2021). Biology. 


Jordan Chapman, Ahmed E. Ismail and Cerasela Zoica Dinu (2018). Industrial Appli- 
cations of Enzymes: Recent Advances, Techniques, and Outlooks; catalysts. 


Mishra, Rashmi. (2020). Optimization of culture parameters for a-glucosidase pro- 
duction from suspension culture of moss Hyophilla nymaniana (Fleish.) Menzel. 
Journal of Genetic Engineering and Biotechnology. 18. 


Malhotra, G., Chapadgaonkar, S.S. (2020). Taguchi optimization and scale up of 
xylanase from Bacillus licheniformis isolated from hot water geyser. J Genet Eng 
Biotechnol 18,65. 


Mishra, R., Chandra, R. (2020) Optimization of culture parameters for a-glucosidase 
production from suspension culture of moss Hyophilla nymaniana (Fleish.) Menzel. 
J Genet Eng Biotechnol 18, 82. 


Hammami, A., Bayoudh, A., Abdelhedi, O. (2018). Low-cost culture medium for the 
production of proteases by Bacillus mojavensis SA and their potential use for the prep- 


aration of antioxidant protein hydrolysate from meat sausage by-products. Ann 
Microbiol 68, 473-484. 


Rouaa Daou, Karine Joubrane, Richard G. Maroun, Lydia Rabbaa Khabbaz, Ali 
Ismail, André El Khoury (2021). Mycotoxins: Factors influencing production and 
control strategies, AIMS Agriculture and Food; Volume 6, Issue 1: 416-447 


Abiola Ezekiel Taiwo, Tafirenyika Nyamayaro Madzimbamuto, Tunde Victor 
Ojumu, (2020). Optimization of process variables for acetoin production in a 
bioreactor using Taguchi orthogonal array design, Heliyon, Volume 6, Issue 10. 


Sunil Chamoli, (2015). ANN and RSM approach for modeling and optimization of 
designing parameters for a V down perforated baffle roughened rectangular channel, 
Alexandria Engineering Journal, Volume 54, Issue 3, Pages 429-446. 


International Education & Research Journal [IERJ] 


E-ISSN No : 2454-9916 | Volume: 8 | Issue: 6 | Jun 2022 


